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1. Introduction 

The interaction between the gravitational field and an elastic solid body became quite 
an important problem in astrophysical applications, since the discovery that the crust of 
neutron stars probably exists in the form of a solid, due to the process of crystallization of 
dense neutron matter (see e.g. McDermott, Van Horn & Hansen 1988). Recently, the present 
authors proposed a new approach to relativistic elasticity (Kijowski & Magli 1992). The aim 
of the present paper is to show the relationship between our approach and the theory proposed 
by Carter and Quintana (1972) and developed in Carter (1980). It turns out that the method 
of deriving the dynamical equations of the theory used by the latter authors is analogous to 
the derivation of (special-relativistic) Maxwell equations via the (general-relativistic) Hilbert 
variational principle, without introducing the notion of electromagnetic potential. On the 
other hand, our approach is analogous to the standard description of electrodynamics in 
terms of potentials. 

In our formulation all the physical quantities (like e.g. the stress and the strain tensors, 
the matter current and so on) are defined in terms of first order derivatives of the potentials. 
This way, all the compatibility conditions of the theory are automatically satisfied. As a 
consequence, the dynamics can be formulated in terms of three (independent) second-order 
hyperbolic partial differential equations imposed on three (independent) unknown functions: 
the gauge potentials. This simplifies considerably the dynamical structure of the theory. 

The equivalence between the two formulations is a straightforward consequence of the 
following, remarkable feature of the relativistic mechanics of continuous media: the symmetric 
and the canonical energy-momentum tensors of the theory do coincide (due to the convention 
which is generally used, they actually coincide up to a sign - see e.g. Jezierski & Kijowski 
1991). It is well known that this is not the case for a general relativistic field theory, like e.g. 
electromagnetism, although the relationship between the two tensors is well understood (see 
e.g. Kijowski & Tulczyjew 1979). In the particular 

case of relativistic elasticity, the simplest way to understand the equivalence between 
the two tensors is based on the following observation. Relativistic elasticity interacting with 
the gravitational field may be regarded as the theory of two symmetric tensor fields: the 
physical metric g and the material metric h (see section 2). As will be seen in the sequel, the 
Lagrangian of the theory depends on both metrics via their combination g^'^hj^x only. Hence, 
its variations with respect to both metrics coincide (up to a sign). 

2. Kinematics 

Let M. be the general-relativistic space-time equipped with a pseudo-riemannian metric 
tensor g^j, (^u, z/ = 0, 1, 2, 3) of signature (— , +, +, +). Suppose that M. (or an open domain 
O <Z M.) is filled with a continuous material. The collection of all the trajectories of the 
material points may be considered as a 1-dimensional foliation of O. This foliation can be 
described, as usual, by a normalized, time-like, future oriented vector field w^, tangent to the 
trajectories and called the four-velocity: 

u^u^ = -1 . (1) 

To complete the description of the configuration of the material one has to include the in- 
formation about such quantities as the matter density, internal strains, temperature etc., 
according to the physical character of the material. In the case of a barotropic and isotropic 
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elastic body, which we are going to describe in the present paper, the complete information 
about the configuration is carried by its internal material metric. It describes the "would be" 
rest-frame distances between adjacent molecules of the body, if the corresponding infinites- 
imal portion of the body had been extracted from the rest of the material and left in the 
perfectly relaxed state. Mathematically, the material metric will be described by a symmetric 
tensor field h^i,. 

The "material distance" between points belonging to the same trajectory has to vanish, 
because they correspond to the same molecule of the body. Hence, h^^^ has to be orthogonal 
to the velocity: 

hf^^u^ = . (2) 

This implies that h has the signature (0, -|-, 4-, +), and therefore it may assume the value of 
any non-negative, symmetric tensor, having precisely one vanishing eigenvalue, and such that 
the corresponding eigenvector is time-like. At each point of Ai there is a 9-parameter family 
of such objects. Observe that the information about the velocity u is already contained in 
h. Indeed, u may be defined as the unique time-like, future-oriented, normalized eigenvector 
of h. We conclude that h carries the entire information about the configuration of the 
physical system in question. The dynamical equations of relativistic elasticity theory may be 
formulated as first order differential equations imposed on the quantity h. 

We are going to describe only materials without memory, i.e. such that the material dis- 
tance between two adjacent particles remains constant during the evolution. More precisely, 
this condition means the following: extracting an infinitesimal portion of the material and 
letting it relax leads to the same distance between the particles, independently of the moment 
at which the portion has been extracted. Mathematically, this means that the metric h is 
"frozen" in the material, i.e. its Lie derivative with respect to u vanishes: 

where 

£uhij,u :=u^Vxh^u + hfj^xV^u^ + KxVf,u^ = 

=u^dxh^^ + h^xduU^ + Kxd^u^ . 

As u is also a function of h, equation (3) has to be considered as an identity imposed on h 
alone; only those h describe physically admissible configurations of the material, which fulfill 
this condition. Due to the identity 

there are only 6 independent conditions in (3), imposed on the 9 independent components of 
h. We conclude that the configurations of an elastic material are described by 3 independent 
functions (degrees of freedom) defined implicitly by the identities (3) . 



3. Relativistic strain tensor 

The material is locally relaxed at a point a; G A^ if and only if the physical distances 
described by the metric g coincide with the material distances described by h. This happens 
if the material metric coincides with the physical metric on the subspace orthogonal to u, 
i.e. if the following equation is satisfied at x: 

hfiv = E^i, , (5) 
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where by E we denote the orthogonal projector 

In a generic situation h is not equal to E. The bigger is the difference between them, the 
stronger is the state of strain of the material. There are many ways to measure this state; 
one possibility is to introduce their difference (Cattaneo 1973, Maugin 1978): 

^tiu '■= -^{E^i, — h^i,) . (6) 

Another description has been proposed by the present authors (see Kijowski & Magli 1992) 
in terms of the quantity 

S:=\\ogK (7) 

where 

i^^ := ^'^^/lA. - u^u,) . 

Such strain tensors are both orthogonal to u and vanish for the locally relaxed state. There 
is a one-to-one correspondence between them. Therefore, both descriptions are equivalent. 
FFrom the theoretical point of view, the description (7) is somewhat preferable since S is 
free to assume any value of a symmetric tensor orthogonal to u, whereas E is subject to a 
rather involved matrix inequality [h — E — 2S > 0). In the linearized version of the theory, 
both descriptions obviously coincide. 

The internal elastic energy accumulated in each portion of the body is a function of 
its deformation, described by its state of strain. Hence, the energy depends upon both the 
physical metric g and the material metric h via their combination S (or E). 

4. Carter- Quintana variational principle 

A kinematically admissible h may describe a real physical situation if and only if it 
satisfies also the dynamical equations of the theory, which can be derived from a variational 
principle. For this purpose Carter and Quintana consider a bigger physical system, composed 
of both the elastic body and the gravitational field interacting with it. The total Lagrangian 
density of such a system is: 

1 



where e is the rest-frame energy density of the material and R is the scalar curvature of the 
space-time. Keeping the material configuration h fixed and varying A with respect to the 
gravitational field g we obtain Einstein equations 

G^, = SttT^, (8) 

where T^i, is the symmetric energy-momentum tensor: 

^M--=2^-f^^-- (9) 



J. Kijowski, G. Magli, Relativistic elastomechanics is a gauge-type theory 5 

In this approach the dynamical equations for h arise only as the compatibility conditions 

V^T'^^ = (10) 

of (8) with the Bianchi identities, satisfied by the Einstein tensor G. Among eqs. (10), only 
three are independent (elasticity theory has three degrees of freedom, as we have seen) and 
in fact in the Carter and Quintana paper it is shown that the equation 

ti^V^T^" = 

holds identically. The proof is rather involved and requires the introduction of the concept 
of convected derivative. 

But the fundamental drawback of the Carter-Quintana method of deriving dynamical 
equations is that elasto dynamics and geometro-dynamics can not be separated. Indeed, the 
gravitational field can not be given a priori, as e.g. in special relativity. The theory of "test" 
elastic bodies in the fiat Minkowski space is therefore excluded, because the space-time metric 
has always to be considered as a dynamical variable. Also the non-relativistic theory, which 
we know to be a Lagrangian theory (see e.g. Sommerferld 1950) is automatically excluded. 
Moreover, it is very difficult to impose upon the general theory the existence of symmetries 
in order to handle specific problems, as the equilibrium of axisymmetric solid stars (see e.g. 
Carter 1973, Carter & Quintana 1975, Quintana 1976, Prion 1992). 

To illustrate the Carter-Quintana approach let us consider classical. Maxwell electrody- 
namics. Here, the field configurations are described by the electromagnetic, skew-symmetric 
tensor F^^. In our example F plays a role analogous to that of h in elasto dynamics. The 
first pair of Maxwell equations: 

df^Fp^ + dpF,^ + d^F^p = (11) 

can be considered as the kinematical condition, analogous to (3). We will show that the 
remaining dynamical Maxwell equations can be obtained from the general-relativistic proce- 
dure, analogous to the Carter-Quintana approach. For this purpose, consider the Lagrangian 
of the system composed of both the electromagnetic and the gravitational field: 

Keeping the electromagnetic field F fixed and varying A with respect to the gravitational 
field g we obtain Einstein equations: 

G;; = Stt (p^PFp^ + \F'^PF„p5'i 

Now, Maxwell equations may be obtained as the compatibility conditions for the above sys- 
tem. Indeed, the Bianchi identities imply 



Fp.VpF'^P + F'^PVpFp, + -FP^V^Fp^ = . 



1 
2" 
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Using the skew-symmetry of F and the kinematical equations (11), it is easy to see that the 
last two terms cancel. Hence, for a generic, non singular F, we obtain the second pair of 
Maxwell equations 

V^F^" = . (13) 

Of course, the above approach to electrodynamics, although equivalent to the standard 
one, is very inconvenient if one wants to describe the electromagnetic field in a given space- 
time geometry (e.g. in Minkowski space). 

5. Relativistic elasticity in terms of potentials. 

To obtain the standard variational principle for the Maxwell field (not necessarily coupled 
with the gravitational field), one introduces the electromagnetic potential A^, such that 

F^. = d^A^ - d^A^ . (14) 

In terms of the potential, electromagnetism becomes a gauge theory: 

equations (11) are automatically satisfied and equations (13) become second-order dy- 
namical equations for the potentials. Such equations can be obtained directly from the 
variation of the electromagnetic Lagrangian with respect to the potentials. There is no need 
to couple electrodynamics to gravity, although it is possible. In the latter case the total 
Lagrangian is again equal to (12). 

In the case of elastomechanics, we propose a similar approach. As "potentials" of the 
theory we take the collection of all the idealized "molecules" of the material, organized in 
an abstract 3-dimensional differential manifold Z, called the material space. The space-time 
configuration of the material is completely described if we specify a mapping Q : M. —^ Z . 
To a given point of the physical space-time (i.e. to a given point of the space and a given 
instant of time) the mapping assigns the "molecule" of the material which passes through 
that particular point at that particular time. Given a coordinate system (^") (a = 1,2,3) 
in Z and a coordinate system [x^) in M., the configuration of the material is described by 
three functions ^" = ^""{x^), which play a role similar to that of the potentials A^ = A^{x^) 
in electrodynamics. 

The above potentials may be used for the description of any continuous material (see 
e.g. Kijowski, Pawlik & Tulczyjew 1979 and Kijowski, Smolski & Gornicka 1990 for the 
description of thermo-hydrodynamics). The particular case of a barotropic, elastic material 
requires the use of a metric structure in the material space (Hernandez 1970, Glass & Winicour 
1972, Kijowski & Magli 1992). We assume, therefore, that the space Z is equipped with a 
riemannian (positive) metric 7a6, the material metric. This metric is frozen in the material 
and is not a dynamical object of the theory. It is given a priori for a given material, like e.g. 
the equation of state for a fiuid. To understand its physical meaning consider an infinitesimal 
portion of the material. This portion will tend spontaneously to a relaxed state when no 
external forces act on it. The metric 7 may now be defined as describing the rest-frame space 
distances between adjacent "molecules" , measured in such a locally relaxed state. There are, 
however, materials possessing no globally relaxed state. This happens if Z is not isometric 
to any hypersurface of M., e.g. when Z is curved and M. is fiat. We see that materials with 
internal stresses can also be described in this way. 

Given the space-time configuration, we introduce the space-time version h of the material 
metric as the geometric pull-back of the metric 7 from Z to Ai, i.e. 

h := (rf0*7 • 
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In terms of coordinates we have: 

V ■■= lahi^tu , (15) 

where we denote ^"^ := d^^"'. The 3x4 matrix (9^^") is caUed the relativistic deformation 
gradient. 

Hence, similarly to electrodynamics, the physical field has been expressed in terms of the 
derivatives of the potentials: formula (15) in elastomechanics is analogous to formula (14) 
in electrodynamics. It solves automatically the kinematical condition (3), just as (14) solves 
automatically the kinematical condition (11) in electrodynamics. Indeed, we have: 

£uh^,u = u^dxh^^ + h^xduU^ + K\d^u^ = u^ {dxh^^ - d^h^x - d^Kx) 



■ {d„'^ab)C^i X - labCfiu^ X - labCf^^ \u 



(16) 



To prove that the above expression vanishes identically, let us first observe that the deforma- 
tion gradient is automatically orthogonal to the velocity vector: 

u^C^ = , 

because u is tangent to the trajectories of the material molecules, i.e to the lines given by the 
equation ^" = const. This implies also: 

u\dxJab)=u'Cx^^O^ 

and the other terms in (16) cancel, which ends the proof. 

The theory is, of course, invariant with respect to reparametrizations of the material 

space, which play the role of gauge transformations. Therefore, the fields ^" may be regarded 

as gauge potentials for the "elasticity field" h. They describe the three degrees of freedom of 

the system in an explicit way. 

The group of gauge transformations of the entire theory (elasticity 

interacting with gravity) is therefore the product of the group of space-time diffeomor- 

phisms (which is the gauge group of general relativity) by the group of diffeomorphisms of 

the material space. 



6. Dynamics 

The physical laws describing the mechanical properties of the elastic material can now 
be formulated in terms of a system of second order hyperbolic partial differential equations 
for the 3 unknown fields ^". The equations can be obtained from the Lagrangian 

A = -v^e . (17) 

considered as a function of the fields (^") and their first derivatives. They assume the form 
of the Euler-Lagrange equations: 
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We are going to prove the equivalence of the above equations with the equations (10) of 
Carter-Quintana. For this purpose, we consider the canonical energy-momentum tensor of 
our theory: 

ri' ■■= -^ f ^r. - 5^a) • (19) 

Due to the standard Nother argument, the Euler-Lagrange field equations imply the energy- 
momentum conservation: 

V^T^'^ = . (20) 

The same argument shows (see Kijowski & Magli, 1992) that the identity 

u^V^T^^ = 

holds automatically because of the gauge invariance of the theory. Hence, in the case of 
elastomechanics there are only 3 independent equations among the 4 equations (20). To 
prove their equivalence with the 3 Euler-Lagrange equations (18), it is now sufficient to show 
that both energy-momentum tensors coincide up to a sign: 

In fact, this is a particular case of the general Belinfante-Rosenfeld theorem (Belinfante 
1940, Rosenfeld 1940). The theorem uses the invariance of the Lagrangian with respect to 
the space-time diffeomorphisms. This means that it may be easily proved in framework 
which is more general than the one used so far in the present paper, namely, we are free to 
consider anisotropic elastic bodies as well. In order to introduce the variational principle for 
such bodies, we define the following tensor in the material space: 



from the geometrical point of view, such a tensor is the "image" of the physical metric in the 
material space, i.e. G"^ plays in Z the same role played by h^^, in M.. For a general, possibly 
anisotropic body, we assume the energy density to be a function of the whole G"'': 

e-e(G'^^e). 

The particular case of isotropic bodies corresponds to an internal energy which is a function 
of the invariants of G"'' only; of course, the invariants of G"'' coincide with the corresponding 
invariants of hfj,u-, like e.g.: 

G^ = 7a.G"^ = ia,e,e.g^^ = Kug^^ = h^, . 

According to the definition (19), we have: 

de _ „„ de dO^'^ _ „„ de 



On the other hand, due to definition (9), we have 

Of „ Of dO"^ „ Of 



'^^'^ ~ ^ dqf^" ^^'^'^ ~ ^ dG"-^ ^q^^'' ^^^^'^ ~ ^ QQab^^'f^^ '^ (^Qf^u , (22) 
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which ends the proof of the Behnfante-Rosenfeld identity (21). It is worthwhile to note that 
the first term in (22) is automatically orthogonal to the velocity. Thus, the energy-momentum 
tensor assumes the canonical form 

where the pressure tensor p^i, is orthogonal to the velocity and is given by 

de 



7. Concluding remarks 

In principle, the theory presented here follows the lines of the approach to relativistic 
elasticity proposed by Cattaneo (1973). There is, however, a fundamental difference: Catta- 
neo describes the 

configuration of the material always with respect to a fixed, static reference configuration. 
Choosing a space-like surface {t = const} of this reference configuration we obtain a specific 
representation of our material space Z. We find it, however, more natural to use only the 
abstract, metric structure of Z, instead of working always with two different space-times: one 
to describe the actual configuration and another to describe the reference configuration. Both 
theories are equivalent in the particular case of fiat materials (no internal-frozen stresses, fiat 
material metric) and fixed gravitational field (e.g. special relativity). In order to study the 
interaction between gravity and elasticity in Cattaneo's formulation, it is necessary to follow 
an ad hoc approach, which consists in defining a reference state of the material when the 
gravity is hypothetically "switched off" (Newton constant set equal to zero), and then in 
studying the evolution during a process of adiabatic restore of the coupling constant to its 
own value (Cattaneo 1973, Cattaneo & Gerardi 1975). Mathematically, this is a nice way 
to overcome the difficulties of defining a reference state. From the physical point of view, 
however, switching off the gravity is poorly justified. 

An approach similar to ours has also been given by Maugin (1978). However, this author 
works in a "direct" picture rather than in the "inverse" picture, proposed in our paper. The 
difference between these two pictures consists in inverting the role of the space parameters 
(x^) and the fields (^"). In the "direct" picture the configuration of the material is described 
by the 3 fields x^ depending upon 4 independent parameters (x°,^"). For this purpose a 
(3-|-l)-decomposition of the space-time has to be chosen, which breaks the explicit relativistic 
invariance of the theory. However, the theory proposed by Maugin remains relativistic- 
covariant and its results are independent of a specific choice of the (3+l)-decomposition. 

The advantage of our approach consists in the fact that we really eliminate all the 
constraints in a fully relativistic-invariant way, reducing the number of independent degrees 
of freedom to the three independent functions ^". We hope this approach to be useful for 
specific applications like e.g. the description of the equilibrium and the oscillations of neutron 
stars crusts in a general relativistic context. Some steps in this direction have already been 
done (Magli & Kijowski 1992, Magli 1993a,b). Moreover, as will be shown in a forthcoming 
paper (Kijowski & Magli 1993), our formulation naturally leads to the hamiltonian version 
of the theory, the canonical variables being the fields ^" and their conjugate momenta tt^ = 
dA/d^"'. The Poisson bracket between them assumes its canonical, delta-like form. The 
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above hamiltonian structure is common for any (relativistic or non-relativistic) mechanics of 
continuous media (see also Kijowski, Smolski & Gornicka 1990 and Jezierski & Kijowski 1991). 
It is given by a standard Legendre transformation, once the Lagrangian (17) is expressed in 
terms of the (independent) variables (^") and their first (unconstrained) derivatives (see also 
Kijowski & Tulczyjew 1979). 
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